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Vertical tail buffeting is a serious multidisciplinary problem that limits the performance of twin -tail fighter aircraft. 
The buffet problem occurs at high angles of attack when the vortical flow breaks down ahead of the vertical tails 
resulting in unsteady and unbalanced pressure loads on the vertical tails. This paper describes a multidisciplinary 
computational investigation for buffet load alleviation of full F/A-18 aircraft using distributed piezoelectric 
actuators. The inboard and outboard surfaces of the vertical tail are equipped with piezoelectric actuators to control 
the buffet responses in the first bending and torsion modes. The electrodynamics of the smart structure are expressed 
with a three-dimensional finite-element model. A single -input-single -output controller is designed to drive the active 
piezoelectric actuators. High-fidelity multidisciplinary analysis modules for the fluid dynamics, structure dynamics, 
electrodynamics of the piezoelectric actuators, fluid -structure interfacing, and grid motion are integrated into a 
multidisciplinary computing environment that controls the temporal synchronization of the analysis modules. Peak 
values of the power spectral density of tail tip acceleration are reduced by as much as 22% in the first bending mode 
and by as much as 82% in the first torsion mode. RMS values of tip acceleration are reduced by as much as 12%. 


NOMENCLATURE 

du longitudinal induced strain coefficients 

dy lateral induced strain coefficients 

Ei longitudinal Young’s modulus (fiber axis) 

E 2 lateral Young’s modulus (electrode axis) 

Ek electric field components 

[e] piezoelectric matrix in stress-charge form 

G modulus of shear stress 

K coefficient of relative permittivity 

K uu structural stiffness matrix 

K u(p coupled electric -elastic stiffness matrix 

Kqy stiffness matrix of the piezoelectric material 

M uu structural mass matrix 

u structural nodal displacement 

v Poisson’s ratio 

(p electrostatic potential 

[e] dielectric matrix 


INTRODUCTION 

T he capability of fighter aircraft to maneuver at 
high angles of attack suffers significantly from the 
multidisciplinary problem of vertical tail buffeting. At 
high angles of attack, the vortices that emanate from the 
leading-edge extension (LEX) of the wing break down 
upstream of the vertical tails. In these conditions, the 
breakdown vortical flow impinges upon the vertical tails 
resulting in severe dynamic buffet loads. 1 * ' 5 The 
breakdown flow contains significant amount of energy 
over a frequency bandwidth that often coincide with the 
predominant structural modes of the vertical tails. The 
buffet loads imposed upon the vertical tails resulted in a 
premature failure of the vertical tails and increased the 
inspection and maintenance cost of fighter aircraft. 

There are essentially two major techniques to control 
the tail buffet problem: flow control or structure control. 
The flow control methods aim into modifying the 
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vortical flowfield around the vertical tails to reduce the 
buffet loads. Passive flow control 6 ' 8 and active flow 
control 9,10 methods have been proposed, but only the 
passive methods received the most attention. As a near 
term to the buffet problem of the F/A-18 aircraft, 
passive stream wise LEX fence devices were used. The 
LEX fences were fitted over the upper surface of the 
wing LEX to reconstruct the vortical flow over the 
aircraft. 6 The fences, although developed through trial 
and error wind-tunnel experiments, reduces the buffet 
loads experienced by the vertical tails. In flight tests by 
Lee et al, 7 the acceleration peaks at a point close to the 
tip of the vertical tail were reduced from 450 to 200 g 
with the addition of the LEX fence. However, a loss of 
3 percent of maximum lift was monitored. The fences 
were also less effective at high angles of attack. Other 
passive flow control methods, such as dorsal fin 
extension of the vertical tails or up-deflected LEX, have 
been proposed by Rao et al. 8 Although these methods 
provided some buffet relief, an aerodynamic penalty 
was imposed on the aircraft. These methods are also not 
effective at all flight conditions. 

Active control of the tail structure can be achieved using 
either active control of the tail rudder 11,12 or using 
piezoelectric actuators. 13,14 In the active rudder method, 
the rudder is deflected in response to a feed back signal 
from the motion of the tail tip. Moses 12 applied the 
active rudder concept to l/6 th F/A-18 model in the TDT 
wind tunnel at NASA Langley. The active rudder was 
successful in reducing the buffet responses at the first 
bending mode around 15 Hz. However, the actuation 
bandwidth of the rudder is limited. For example, the 
rudder bandwidth of the F/A-18 aircraft is less than 20 
Hz. Therefore, it is not effective for reducing the buffet 
responses at the second mode around 45 Hz (first 
torsion mode). 

The piezoelectricity (Greek for pressure electricity) is a 
property of materials that generate an electrical charge 
when the materials are subject to material strain, and 
conversely, piezoelectric materials deform when they 
are subject to an electrical field. A piezoceramic is 
therefore capable of acting as either a sensor or actuator. 
Piezoceramics can be used for control of vibrations by 
staking them where excessive vibrations occur. The 
strain of the structure is fed into a controller that injects 
voltage into the piezoelectric (PZT) patches. The 
mechanical response of the PZT patches due to the input 
voltage injects vibrations into the structure, which is 
equivalent and opposite to the initially detected one 
such that the net vibration becomes minimum. 

As part of the Actively Controlled Response of Buffet- 
Affected Tails (ACROBAT) program, Moses 12,13 
conducted wind tunnel experiments on l/6 th F/A-18 


model. The starboard vertical tail was equipped with 
active rudder, and the port tail was equipped with 
piezoelectric actuators. The power spectral density 
(PSD) of the root bending moment at the frequency of 
the first bending mode was reduced by up to 60%. 

Recently, Moses et al 14 tested new sets of actuators, 
active fiber composite actuators (AFC) and macro fiber 
composite actuators (MFC), for buffet alleviation of 
1 /6-scale model of F/A-18 aircraft. A total of ten 
actuators, five per side, were embedded beneath the 
fiberglass shell of the vertical tail, as shown in Figure 1. 
The new actuators are poled electrically in the plane of 
the actuators rather than through the thickness. The new 
actuators performed superbly in reducing structural 
responses caused by buffet. Additional experiments are 
currently underway to test the effectiveness of the new 
actuators in reducing structural responses of the full- 
scale aircraft. 



Figure 1. F/A-18 vertical tail with PZT actuators. 14 


In this paper, a unique high-fidelity multidisciplinary 
computational investigation of active control of the 
vertical tail buffeting of full F/A-18 aircraft using 
piezoelectric actuators is presented. Several control laws 
and aeroservoelastic models have been tested and their 
representative results are discussed. 


THE TECHNICAL APPROACH 

Active control of tail buffeting is a complex 
multidisciplinary aeroelastic problem that involves 
interaction between several physical and numerical 
disciplines. The physical disciplines include the fluid 
dynamics of the vortical flow field, structure dynamics 
of the flexible bodies, and electrodynamics of the PZT 
actuators. The numerical disciplines include grid 
deformation and fluid -structure interface coupling. In 
this investigation, the computational modules of these 
disciplines are integrated into the multidisciplinary 
computing environment (MDICE). 15 MDICE controls 


2 

American Institute of Aeronautics and Astronautics 



the temporal synchronization of the data transfer 
between the multidisciplinary modules. The particular 
sets of analysis modules used in the current 
investigation are outlined next. 

The Fluid Dynamics Module 

The buffet problem of vertical tails occurs at high 
angles of attack where the flow field is characterized by 
strong vortical flow and massive three-dimensional flow 
separation. Therefore, the problem mandates the 
solution of the full Navier- Stokes equations to account 
for vorticity evolution, convection, shedding, and vortex 
breakdown. The full Navier-Stokes equations are solved 
using CFD-FASTRAN 16 module using an implicit 
finite-volume upwind scheme with Roe's flux-difference 
splitting for spatial differencing and a fully implicit 
upwind scheme for temporal differencing. On the 
aircraft surfaces, no-slip and no-penetration conditions 
are enforced. Riemann's inflow/outflow boundary 
conditions are assumed at the far field boundaries in all 
directions. The flow field around the aircraft is assumed 
symmetric, and therefore, a plane of symmetric 
boundary condition is used at the symmetry plane. 

The Structure Dynamics Module 

The structural dynamics equations of the flexible 
vertical tail under the influence of the aerodynamics 
loads of the vortex breakdown flowfield are given by 

[mJu}+[c u Ju}+\kJ{u}={f] ( 1 ) 

where {u} is the displacement vector, [M uu \ is the mass 
matrix, [C uu \ is the damping matrix, [K uu ] is the stiffness 
matrix, and {F} is the aerodynamics force vector. 


{5} = [5 £ ]{r}+[^f{£} (2) 

{D} = [d]{T} + [e T ]{E}, (3) 

where [ s E ] is the compliance matrix, [d] is piezoelectric 
coupling matrix, £ 7 ] is the dielectric constant matrix 
(permittivity components), {S} is the strain vector, {T} 
is the stress vector, {£>} is the charge density vector, and 
{E} is the electric field vector. Equation (2) is the 
sensor equation, which relates the charge developed in 
the piezoelectric material to stresses or strains. Equation 
(3) is the actuator equation, which relates the applied 
electric field to the strain developed in the material. 


The electric field components, Ek, are related to the 
electrostatic potential / by the relations: 

E k=- V k<P' * = 1 ’ 2 ’ 3 * W 


When the variational principle is applied to both the 
piezoelectric equations and the equations of the solid 
structure, the linear dynamic equations of the coupled 
system reduce to the following system of equations: 
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where u, $ F, and G are the global nodal displacement, 
electric potential, force, and applied charge vectors, 
respectively. K m represents the stiffness matrix of the 
original structure. K u( p represents the coupled stiffness 
between electrical field and elastic strain, and K <$ 
represents the stiffness matrix of the piezoelectric 
material. These stiffness matrices are given as follows: 


The structural dynamics equations of the vertical tail are 
solved using three-dimensional direct finite-element 
analysis. The vertical tail of the F/A-18 aircraft is 
modeled using second order hexahedral elements. The 
vertical tail is the only flexible structure in the aircraft. 
The initial conditions are the undeformed vertical tail. 
The vertical tail is assumed clamped at the root. No 
structural damping is assumed in this investigation. The 
material of the vertical tail is assumed to be Aluminum 
and isotropic. The Young's modulus of elasticity is 
4.Cfc><10 10 N/rr?, the density is 2765 kg/ir?, and the 
Poisson's ratio is 0.33. The first structural mode of the 
vertical tail is first bending mode at 15.5 Hz. The 
second mode is first torsion mode at 46 Hz, and the 
third mode is the second bending mode at 54 Hz. 

Finite -Element modeling of PZT Actuators 

Piezoelectric materials combine the constitutive 
equations of solids and dielectric materials into one 
coupled equation. The linear piezoelectric constitutive 
equations (in strain -charge form) are given by 17 


[k U <p\i= J BleB (p dV 
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where N u and Afy are the shape function matrices for the 
displacement and electric fields. 


If augmented entities are defined for the global 
displacement field, mass matrix, stiffness matrix, and 
applied forces as follows: 



( 11 ) 


then, equation (5) may be written as 


Mu + Ku = F, 


( 12 ) 


Equation (12) is similar to the standard structural 
dynamics equations of the tail structure, equation (1). 
Therefore, the modifications of the solid finite-element 
module become as follows: (1) Introduce additional 
nodal dependent variable $ and (2) Modify mass and 
stiffness matrix accordingly. 


The Fluid-Structure Interfacing Module 

The fluid -structure interfacing algorithm is used to 
project the forces from the fluid flow to equivalent 
forces and moments on the flexible -body structure and 
to feed back the aeroelastic deflections of the structure 
nodes to the flow field grid cells. A conservative and 
consistent interfacing algorithm is used in this 
investigation. Conservative interfaces aim to conserve 
the forces and moments in the interpolation process 
between two grids. Consistency, or virtual work 
conservation, requires that the virtual work performed 
by the solid interface is equivalent to the virtual work 
performed by the fluid interface. The interfacing is 
formulated in the most general sense for maximum 
flexibility. There are no inherent assumptions that the 
fluid grid is matched with the structure grid, either 
through different mesh densities or mesh architecture. 


The Grid Motion Module 

The computational grid is deformed every fluid - 
structure data transfer to accommodate the deformed 
shape of the tail. The six outer boundary surfaces of the 
deformed grid block are kept fixed. The grid is 
deformed using a TransFinite Interpolation algorithm 
(TFI). 18 TFI is an interpolation procedure that deforms 
grids conforming to specified boundaries and it is very 
computationally efficient. The spacing between grid 
points is controlled by blending functions that specify 
how far into the original grid the effect of the new 
position of the flexible body surfaces is carried. The TFI 
routine is invoked automatically when a fluid -structure 
interface is exchanged between application modules. 


The Multidisciplinary Computing Environment 

The multidisciplinary modules used in the current 
investigation are integrated into the multidisciplinary 
computing environment MDICE. 15 MDICE is an object- 
oriented environment that enables the analysis modules 
to run concurrently and cooperatively on a distributed 
network of computers. Using MDICE environment, one 
can avoid giant monolithic codes that attempt to provide 
all modules in a single large computer program. Such 
large programs are difficult to develop and maintain and 
by their nature cannot contain up-to-date technology. 


COMPUTATIONAL AEROSERVOELASTIC 
MODEL OF F/A-18 AIRCRAFT 

The computational aeroservoelastic model of the F/A-18 
aircraft consists of three main parts. The first part is the 
geometry modeling of a full F/A-18 aircraft. The second 
part is a piezoelectric aeroservoelastic model of the 
vertical tail. The third part is a control law to drive the 
piezoelectric actuators. 

F/A-18 Geometry Modeling 

The geometry of the full F/A-18 aircraft is modeled 
using a structured body-fitted grid system with one-to- 
one grid connectivity. The surface grid of the F/A-18 
aircraft is shown in Figure 2. The grid is a multiblock 
H-H grid consisting of 53 structured blocks utilizing 
2.56 million cells (2.68 million grid points). The global 
grid extends 5 wing chords upstream and downstream of 
the aircraft, and 8 wing chords normal to the aircraft. 
The grid is clustered near the apex and leading edge of 
the wing LEX in order to produce robust and well- 
defined leading-edge vortices. 



Figure 2. Surface grid of the F/A-18 aircraft. 


For the structure dynamics analysis, the vertical tail of 
the F/A-18 aircraft is modeled using second -order brick 
(hexahedral) elements. The computational structure 
dynamics (CSD) and computational fluid dynamics 
(CFD) surface grids of the vertical tail are shown in 
Figure 3. The CFD and CSD surface grids of the 
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vertical tail are not identical. A total of 1840 grid cells 
define the CFD surface grid of the tail, while 576 grid 
cells define the CSD surface grid of the tail. The forces 
and deflections are interpolated between the CFD and 
CSD grids using the interfacing module outlined before. 



Figure 3. Surface grid of the vertical tail. 


The PZT Aeroservoelastic Model 

The inboard and outboard surfaces of the 
aeroservoelastic model of the F/A-18 vertical tail are 
shown in Figure 4. A total of five EZT actuators are 
shown on every side of the vertical tail. The root PZT 
actuators (3 per side) were strategically located near the 
tail root to alleviate the buffeting in the first bending 
mode. The tip PZT actuators (2 per side) were located 
near the tail tip to alleviate the buffeting in the first 
torsion mode. The two stacks of actuators on both sides 
of the tail are set in an opposite phase to each other to 
generate bending and torsion. In other words, when one 
side of the tail is contracting the other side of the tail is 
expanding to generate either bending or torsion. Each 
set of PZT actuators on both sides of the tail is being 
commanded by a control law to strain in opposite 
directions simultaneously to enable the vertical tail to 
bend. The acceleration at a near-tip point (75% span and 
80% chord) is recorded and used as the input signal of 
the control law. 


(a) Inboard surface (b) Outboard surface 

Figure 4. Aeroservoelastic model of the vertical tail 

The ceramic PZT patches used in this computational 
study are assumed to be the Macro Fiber Composites 
(MFC). 19 The PZT fibers (PZT- 5 A) of the MFC are 
transversely aligned in order to sense and actuate in- 
plane stresses and strains for structural actuation. The 
material properties of the MFC were obtained from 
Wilkie et al 19 and Azzouz et al. 20 The material 
properties are given in Table 1. 


Table 1. Material properties of the PZT patches. 


Properties 

MFC 

units 

dll/d3i 

530 E- 12 

m/V 

dl2/d32 

-210 E-12 

m/V 

El 

3.65 E+10 

N/m? 

E2 

0.76 E+10 

N/m 2 

G12 

1.46 E+10 

N/m 2 

Q23 

0.73 E+10 

N/m 2 

Kll 

495 

— 

V12 

0.25 

— 

P 

7552 

Kg/rn^ 

Thickness 

0.13 

inches 


The Control Law 

The block diagram of the active control system used for 
tail buffeting alleviation is shown in Figure 5. The 
control law is a single -input/single -output control law 
that utilizes a feedback signal from an accelerometer 
mounted near the tip of the vertical tail. 

The measured acceleration of the vertical tail, u , is fed 
through the controller into two control laws, (CLaw) t i P 
and (CLaw>oot of the tip and root PZT patches, 
respectively. Each of the control laws has different 
control law gain, which is used to set the peak 
magnitude of the control law to unity. A switch is used 
to close or open the control loop of either the tip or the 
root control law. The outputs of the controllers, \i P and 
I root , are sent to a summing junction where it may be 
combined with other signals for system identification. 
The signal is then sent into voltage amplifiers with gains 
of Kip and Koot- The voltage amplifiers are used to 
increase the magnitude of the voltage signal to a value 
within the operational range of the actuators. The 
voltage signal is then applied to drive the PZT actuators 
causing strain actuation. The input to the vertical tail is 
therefore the strain actuation in addition to the 
aerodynamics buffet loads. 

The aerodynamic buffet loads are being input to the 
vertical tails through the unsteady surface pressure on 
both sides of the vertical tail. The root mean square 
(RMS) of the differential pressure on the surface of the 
vertical tail at a point located at 45% chord and 60% 
span is shown in Figure 6. The results are compared 
with flight and wind tunnel data reported by Meyn et 
al. 21 As shown in the figure, the maximum buffet 
condition occurs around the 30° angle of attack. The 
power spectral density (PSD) of the differential 
pressure at the 30° angle of attack is also shown in 
Figure 6. The peak value of the PSD occurs at 16 Hz, 
which is very close to the first bending mode of the 
vertical tail structure . 
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Figure 5 . Active control system diagram for tail buffet alleviation. 


i 

a 




i 

u 


0.4 , 


0.3 


0.2 


O Flight Data (Mtyn vl ii„ 1 W4)| 
q @ D □lull Scale (Meyn ct ttl* 1994) 

n Q ^ Scale (Meyn et ul., 1994) 

• Present Computat km 


B 


B 




(P 




20 


30 40 50 

—Angl e o f Atta ck 


60 


n 



Figure 6. The RMS and PSD of differential pressure at 45% chord and 60% span point of the vertical tail. 



(a) CLaw root ^br the root PZT actuators (b) CLaw lip /or the tip PZT actuators 

Figure 7. Frequency response of the controllers , CLaw roo t and CLa\\fc p . 
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The frequency responses of the control laws, CLawfoot 
and CLaw t jp, are shown in Figure 7. The peak 
magnitude is intentionally set to unity so that 
adjustments to the magnitude could be made via the 
gain block K,. 00l and Kt ip . The phase of the control law 
near the modal frequency of interest is set to -90° to 
augment structural damping. 

Two control laws were designed for the root and tip 
PZT actuators to control tail buffeting in the first 
bending and first torsion modes. The single -input 
single -output controllers were designed using low pass 
filtering for the root controller and band pass filtering 
for the tip controller. The filtering are used so as not to 
excite the other modes when feedback is turned on as 
well as to concentrate the amplifier energy near the 
modal frequencies of interest. 


RESULTS AND DISCUSSION 

The computational experiments of the effectiveness of 
the PZT actuators in alleviation of tail buffeting of full- 
scale F/A-18 aircraft are conducted near maximum 
buffet condition. At Mach number of 0.243 and 
Reynolds number of 11 million, the maximum buffet 
condition for the full-scale F/A-18 aircraft was found to 
occur around 30° angle of attack. 16 Several variations of 
the control laws and different aeroservoelastic models 
were tested computationally. The power spectral density 
(PSD) of the tail tip accelerations is used to compare 
open-loop and closed-loop results . 

This multidisciplinary aeroelastic problem is solved 
using two consecutive steps: the first step solves for the 
vortical flow field characteristics around the rigid 
aircraft configuration. The initial condition of this step 
is the undistributed free-stream condition everywhere. 
The solution is carried out until the changes in the 
vortical flow field become insignificant. The second 
step solves for the unsteady aerodynamic flowfield and 
aeroelastic responses of the flexible configuration. The 
initial condition of this step is the final solution of the 
first step. In the second step, strong coupling between 
the fluid and structure is considered. Strong coupling 
occurs when the inertial effect of the motion of the solid 
is fed back into the fluid. The time steps of the 
computations are fixed in all cases at 1 O' 4 seconds. The 
problem is solved on a Linux computer cluster of six 
units. The total grid is divided almost evenly between 
the six processors using domain decomposition. The 
speed of each processor is 1.0 GHz. Every partition of 
the grid requires approximately 160 MB of memory. 
The computational time is about 40 psec per time step 
per grid point. 


Buffeting Alleviation Results 

In this part of the investigation, the aeroservoelastic 
model shown in Figure 4 is used along with the control 
law shown in Figure 7 to control the buffet problem. 

The effects of closed-loop control (feedback on), 
applied only to the root PZT actuators, on the power 
spectral density (PSD) of tail tip accelerations are 
shown in Figure 8. A maximum input voltage command 
of 2000 volts is applied to the root PZT actuators. The 
peak magnitude of the closed-loop (feedback on) 
buffeting is lower than the peak magnitude of the open- 
loop (feedback off) buffeting for frequencies near the 
first bending and torsion modes. The PSD peak 
magnitude of the tip acceleration in the first bending 
mode is reduced by about 18%. The RMS value of tip 
acceleration is reduced by about 4%. The PSD peak 
magnitude of the tip acceleration in the first torsion 
mode is also reduced by about 9.5 %, although no input 
command is applied to the tip PZT actuators. This 
reduction is caused by excitation at the first torsion 
mode caused by the feedback controller at the first 
bending mode. The increase in magnitude around 47.5 
Hz is caused by the control law and could have been 
prevented by modifying the control law to adjust the 
phase in that region of the spectra. 

The PSD of the tip acceleration is shown in Figure 9 
when closed-loop control (feedback on) is applied to 
both the root and tip PZT actuators. Maximum input 
voltage commands of 2000 volts are applied to both the 
root and tip PZT actuators. The peak magnitude of the 
closed-loop buffeting in the first torsion mode is 
considerably lower than the peak magnitude of the 
open-loop buffeting. The PSD peak magnitude of the tip 
acceleration in the first torsion mode is reduced by 
about 69%. The RMS value of tip acceleration is 
reduced by about 1 0%. The PSD peak magnitude of the 
tip acceleration in the first bending mode is also reduced 
by about 22%. This clearly shows that the PZT actuators 
are much more effective in reducing the peak 
magnitudes at the first torsion mode than those at the 
first bending mode. 

The reason that the PZT actuators are less effective for 
the first bending mode can be explained from the 
distribution of the strain energy and the stiffness of the 
tail structure. In the first bending mode, the strain 
energy of the vertical tail is concentrated near the root 
of the tail. The PZT actuators are not as effective near 
the tail root because of the large structural stiffness near 
the root. In the first torsion mode, the strain energy is 
concentrated near the upper third of the tail where the 
structural stiffness is significantly lower than the 
stiffness at the tail root. 
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Deflection (m) 



(a) PSD of tip acceleration at bending mode (b) PSD of tip acceleration at torsion mode 

Figure 8. PSD of tip acceleration using feedback control to the root PZT actuators only. 



(a) PSD of tip acceleration at bending mode (b) PSD of tip acceleration at torsion mode 

Figure 9. PSD of tip acceleration using feedback control to the root and tip PZT actuators. 



Figure 10. The time history of the tip deflection Figure 11. The time histories of tip acceleration and 

using feedback control to the root and tip PZT the input signal to the PZT actuators. 

actuators. 
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The time histories of the tip deflection at the 75% 
span and 80% chord point of the vertical tail are 
shown in Figure 10 for both the open loop and closed 
loop cases. The closed loop controller is shown to 
reduce the amplitude of the tip deflection. The figure 
also shows that it takes almost half a second before 
the effect of the closed loop controller shows up. The 
time history of the tip acceleration (input signal to the 
control laws) and the time histories of the input 
command signal to the root and tip PZT actuators 
(output signal of the control laws) in volts are shown 
in Figure 11. A maximum of almost 2000 volts are 
kept for both signals. 

Effectiveness of Integrator Control Law 

In this part of the investigation, the effectiveness of a 
simple integrator control law is discussed. In the 
integrator control law, the transfer function of the 
control law is simply set to G(s) = 1/s. The reason the 
integrator control law is tested here is that a pure 
integrator will give an upper bound to the active 
damping, especially when an accelerometer signal is 
fed back to the controller. The magnitude of an 
integrator drops off linearly with frequency, while 
providing constant -90° phase lag for all frequencies. 
However, depending on the frequencies of the other 
modes, there may be sufficient gain to react 
considerably to those other frequencies. In practice, 
integrator control law would only be used when the 
location of the accelerometer naturally isolates the 
mode of interest from the total multi-mode 
acceleration signal. The reason is that the integrator 
can become saturated quickly if it tries to respond to 
all modes, especially those at high frequencies. 

The PSD of the tip acceleration is shown in Figure 12 
when closed-loop integrator controllers are applied to 
both the root and tip PZT actuators. The PSD peak 
magnitude of the tip acceleration in the first bending 
mode is reduced by about 7%. The RMS value of tip 
acceleration is reduced by about 12%. The PSD peak 
magnitude of the tip acceleration in the first torsion 
mode is also reduced by about 38%. 

Alternative PZT Setup of the Aeroservoelastic 
Model 

The patches of the PZT actuators presented in the 
previous aeroservoelastic model, shown in Figure 4, 
were to resemble the configuration setup of the 1 /6 th 
model of the F/A-18 aircraft tested at NASA 
Langley, shown in Figure 1. However, for the full- 
scale model, the size of the patches may be too large 
to be manufactured. Therefore, an alternative setup of 
the PZT actuators of the aeroservoelastic model of 
the F/A-18 vertical tail was assumed using smaller- 


size PZT patches. In this model, shown in Figure 13, 
96 PZT patches (48 per side) were distributed over 
the inboard and outboard surfaces of the vertical tail. 
The root PZT actuators (36 per side) and the tip PZT 
actuators (12 per side) were directed in the same way 
as the other aeroservoelastic model as to damp 
structural vibrations in the first bending and torsion 
modes, respectively. The PZT patches are poled 
electrically along the plane of the actuators, similar to 
the previous model. The control laws that are used to 
drive these PZT actuators are similar to those used 
with the previous model and represented in Figure 7. 

The PSD of the tip acceleration is shown in Figure 14 
when closed-loop controllers are applied to both the 
root and tip PZT actuators. Maximum input voltage 
commands of 2000 volts are applied to both the root 
and tip PZT actuators, similar to the other 
aeroservoelastic model. The PSD peak magnitude of 
the tip acceleration in the first torsion mode is 
reduced by about 82%. The RMS value of tip 
acceleration is reduced by about 12%. The PSD peak 
magnitude of the tip acceleration in the first bending 
mode is also reduced by about 7%. This shows that 
this PZT setup of the aeroservoelastic model is more 
effective than the other model in reducing the peak 
magnitude of acceleration in the first torsion mode. 
However, it is less effective in the first bending 
mode. One reason for the poor performance at the 
first bending mode is that the PZT patches near the 
root of the vertical tail are located in higher locations 
than those of the previous setup of Figure 4. 


CONCLUSION 

An active smart material control system environment 
has been developed to control vertical tail buffeting 
of fighter aircraft. A computational investigation has 
been conducted to control tail-buffeting responses of 
full-scale F/A-18 aircraft using piezoelectric (PZT) 
actuators at maximum buffet conditions. Several 
aeroservoelastic models and several control laws 
have been computationally tested. The PZT actuators 
were strategically located over the inboard and 
outboard surfaces of the vertical tail to alleviate the 
tail buffeting in the first bending and torsion modes. 
The investigation revealed that the PZT actuators are 
more effective in reducing the structural responses in 
the first torsion mode than those in the first bending 
mode. The PZT actuators reduced the power spectral 
density peaks of the tip acceleration by up to 22% in 
the first bending mode and by up to 82% in the first 
torsion mode. Root mean square values of tip 
acceleration are reduced by as much as 12%. 
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(a) PSD of tip acceleration at bending mode. (b) PSD of tip acceleration at torsion mode. 


Figure 12. PSD of tip acceleration using integrator feedback controller to the root and tip PZT actuators. 



(a) Inboard surface (b) Outboard surface 


Figure 13. An alternative PZT setup of the aeroservoelastic model of the F/A-18 vertical tail. 



(a) PSD of tip acceleration at bending mode (b) PSD of tip acceleration at torsion mode 

Figure 1 4. PSD of tip acceleration using feedback control to the root and tip PZT actuators. 
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